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ABSTRACT 

We report on a program of spectroscopic observations of gravitationally-lensed QSOs with 
multiple images. We seek to establish whether microlensing is occurring in each QSO image us- 
ing only single-epoch observations. We calculate flux ratios for the cores of emission lines in image 
pairs to set a baseline for no microlensing. The offset of the continuum flux ratios relative to this 
baseline yields the microlensing magnification free from extinction, as extinction affects the con- 
tinuum and the lines equally. When we find chromatic microlensing, we attempt to constrain the 
size of the QSO accretion disk. SDSSJ1004+4112 and HE1104-1805 show chromatic microlensing 
with amplitudes 0.2 < |Ato| < 0.6 and 0.2 < |Am| < 0.4 mag, respectively. Modeling the accre- 
tion disk with a Gaussian source (/ oc exp(— R 2 /2r 2 )) of size r s oc X p and using magnification 
maps to simulate microlensing we find r s (A3363) = 7 ± 3 light — days(18.f ± 7.8 x f0 15 cm) and 
p = 1.1 ± 0.4 for SDSS1004+4112, and r s (A3363) = 6 ± 2 light - days(15.5 ± 5.2 x 10 15 cm) and 
p = 0.7±0.1 for HE1104-1805. For SDSSJ1029+2623 we find strong chromaticity of - 0.4 mag in 
the continuum flux ratio, which probably arises from microlensing although not all the available 
data fit within this explanation. For Q0957+561 we measure B-A magnitude differences of 0.4 
mag, much greater than the ~0.05 mag amplitude usually inferred from lightcurve variability. It 
may substantially modify the current interpretations of microlensing in this system, likely favor- 
ing the hypothesis of smaller sources and/or larger microdeflectors. For HS0818+1227, our data 
yield posible evidence of microlensing. 

Subject headings: gravitational lensing: strong - gravitational lensing: micro - accretion disks - quasars: 
individual: HS0818+1227, Q0957+561, SDSS1004+4112, SDSS1029+2623, HE1104-1805 
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1. Introduction 

Gravitational lenses are a powerful tool to 
study not only the structure of the lensed quasar 
but also the composition of the lens galaxy 
(ISchneider et al.ll992 ; Kochanekl2004 ; Wambsganss 



20011 ICollev et al l 120021 ISchechter etldl 12003 ; 



Fohlmeister et al.l 120081 ). leading to several years 



2006). Simple lens models arc usually sufficient 



to reproduce the positions of lensed QSO images, 
but they can fail to reproduce the optical fluxes of 
these images. The so-called flux ratio anomalies 
are thought to be produced by small-scale struc- 
tures in the gravitational potential of lens g alaxies 



of monitoring. 

On the other hand, most of the microlensing 
searches using optical imaging have been concen- 
trated on quadruple lenses because the effect of 
subst ructure is more important at high magnifica- 
tion ( Witt et al. [l995t ISchechter fe Wambsganssl 
120021: IPoolev et al.ll2007t) . Thev provide enough 
constraints to fit the simplest singular isother- 
mal (SIS) model (ISchechter fc Wambsganssl 1 2003 



dWitt et al.l ll995t iMao fc Schneiderl Il998l: IChiba 



Kochanek fc Dalalll2004l ) and find the flux anoma- 



2002t lMetcalf fc Madeaull200lHDalal fc Kochanek 
!3chechter fc Wambsganssl 2002t Keeton 



2002; 



2002[ iBradac et al.l l2002t iMetcalf fc Zad l2002t 

Moustakas fc MetcalJ 2003 ) . These structures are 
either dark matter subhalos or stars and the ef- 
fects they produce are referred to as millilensing 
and microlensing respectively. 

A substructure is able to produce a flux 
anomaly if the radius of its Einstein ring is large 
compared to the emitting region. Since the sizes 
of the quasar continuum emitting regions depend 
on the wavelength, microlensing by stars in a lens 
galaxy will yield a wavel ength-dependent magnifi- 
cation of the continuum (IWambsganss fc Paczvhski 

1991: Wisotz kTet al.ll9 95: Mos auera. Munoz fc Mediavilla 

^iiu^ ' 

20091 : iMediavilla et al.ll201ll ) that can be strong for 
the UV and the optical but is negligible for the IR. 
Microlensing could also affect the high ionization 
broad emission lines (BEL) that are expected to 
arise from the inner part of the broad-line region 
(BLR). Specifically, microlensing would affect the 
broad wings of the profiles of the high ioniza- 
tion lines that correspond to high velocity emit 



lies. Double lenses, however, do not provide 
enough constraints to produce such a model 
unless the fluxes are used as additional con- 
straints. I n these cases, the flux ratio of emis 
sion l ines dWisotzki et al.l 11993c 



Mediavilla et al 



al. 
17) 



I2009L l201lh or in the infrared (|Agol et al.l 1200 
has been used assuming that the emission regions 
are larger than the microlensing source and dust 
extinction is negligible. 

In the present paper, we use spectra of lensed 
quasars as an alternative approach to photometric 
monitoring, to study microlensing in quadruple or 
double lenses. For each pair of images of a lensed 
quasar (A and B), we base our analysis on the 
measurement of the offsets of the flux ratio of the 
core of the emission lines (Fb/ Fa) core compared 
to the flux ratio of the continuum (Fb / FA) C ont- 
This analysis allows us to distinguish between mi- 
crolensing and dust extinction wi thout assuming a 



mode l for the lensed system (see IMediavilla et al 



ters, l eaving unchanged t h e core (IPopovic et al 



2001; Abaias et al. 
Lewis fc Ibatal 12004 : 



2002; Richards et al 
Gomez- Alvarez et al. 



200 



200 



Low ionization BEL and narrow emission lines 
(NEL) arise from considerably larger regions and 
are supposed to b e insensitive to microlensing 
(jAbaias et al.l l2002h although the low ionization 
BEL of some low-luminosity lensed AGNs may be 
slightly affected by microlensing. 

Observational studies aimed at measuring the 
microlensing effect in lensed QSOs usually con- 
sist of broadband observations repeated over ex- 
tended period s, longer than t he time delays for 



image pairs (jWozniak et al. 2000; lOscoz et al 



201 lL and references therein). In this way, a 
single-epoch spectroscopic observation can suffice 
-through the measurement of microlensing- to es- 
timate physical parameters of interest of the lens 
galaxy (like the fraction of mass in compact ob- 
jects, Mediavilla et al. 2009) or of the unresolved 
quasar source (like the size or the radial tempera- 
ture profile). 

In section 2 we present the data for 5 gravitationally- 
lensed quasars with multiple images for which 
we have obtained low-resolution spectra with 
signal-to-noise greater than 40. The systems, 
HS0818+1227, Q0957+561, SDSS1004+4112, 
SDSS1029+2623, and HE1104-1805, were selected 
because the separations between the images were 
larger than 3". Section 3 is devoted to present the 
data analysis methodology. We discuss our results 
in §4 and give some concluding remarks in section 
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§5. 



Observations and data reduction 



Microlensing detection using spectra of lensed 
QSOs has stringent requirements. First of all, 
we need high signal-to-noise ratio in the spec- 
tra (SNR ~ 40) and sufficient spectral resolution 
(v res < 200 km s — ) to resolve the shape of the line 
profiles. Second, we need to obtain simultaneous, 
spatially-separated spectra of pairs of lensed QSO 
images to compare their continuum and emission 
lines at different wavelengths, which requires good 
seeing conditions (< 0.8"). These requirements 
are achieved using 6-8 m class telescopes under 
good seeing conditions. 

We observed the sample on 11 and 12 January 
2008 with the Blue Channel spectrograph on the 
MMT. Table Q] shows the log of observations. We 
also observed HE1104-1805 on 07 April 2008 with 
the FORS2 spectrograph at the Very Large Tele- 
scope (VLT). Our ground-based observations were 
acquired under excellent atmospheric conditions 
(Table Q}. For Q0957+561, we used archival dat£0 
obtained with the STIS spectrograph on the Hub- 
ble Space Telescope (HST). Components A and B 
were observed with HST at different epochs to ac- 
count for time delay variations in the continuum 
spectra (see Table [T]). A detailed description of 
these obser vations and the sp ectrum analysis can 
be found in iHutchingsl (|2003l ). 

We performed the data reduction with irafU 
tasks. These included bias subtraction, flat field- 
ing, extraction of 1-D spectra and wavelength cal- 
ibration. As we are currently interested only in 
flux ratios (i.e. magnitude difference tub — ttla = 
— 2.b\og(FB / Fa)), we did not flux-calibrate our 
data. Cosmic-ray rejection was carried out in 
those cases where we had at least three exposures. 
The data obtained from the HST archive are al- 
ready fully reduced. 

In spite of the careful data reduction process, 
some systematic errors can affect our measure- 



1 archive data were obtained at the Space Telescope Science 
Institute, operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 
5-26555 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation 



ments. We discuss these in the following para- 
graphs. 

2.1. Spectrum Cross-Contamination 

To avoid cross-contamination between the spec- 
tra of lensed QSO image pairs we selected pairs 
with separation much wider than the average see- 
ing (< Of! 7), which was larger in turn than our 
typical seeing. Our pairs have separations ranging 
from 2'.'6 to 22'.'6 as shown in Table [TJ we esti- 
mate that cross-contamination is negligible in our 
observations. 

2.2. Long- Slit Losses 

To obtain simultaneous pairs of spectra with a 
single slit, we did not observe our targets at the 
parallactic angle; in each case we used the po- 
sition angle defined by the two components. In 
these ground-based observations, we lost a small 
amount of the blue part of the spectra because 
part of the blue quasar light may fall outside the 
slit. Our airmass range was 1.02 — 1.77; we used a 
l'.'O slit width. We used a program developed by 
E. Marchetti at ESCd to calculate the differential 
atmospheric refraction (DAR) for each given wave- 
length and airmass. The atmospheric parameters 
(temperature, humidity, and pressure) we used are 
T = 11.5 °C, H = 14.5 %, P = 743 mbar for VLT 
and T = 12.0 °C, H = 13.4 %, P = 741.6 mbar 
for MMT. The focus for FORS2 is set at 5000 A 
which is our reference wavelength to calculate the 
displacement due to DAR. We estimated the rela- 
tive loss in the flux of each pair of images at 3500 
A and 8000 A, and at 4500 A and 9500 A for 
the MMT and VLT data respectively. We found 
relative losses of < 1% in all MMT spectra; for 
the VLT spectra, we found losses of < 9% due to 
an error in the position angle used in the observa- 
tions. Considering the separation between pairs of 
images (2'.'6 to 22'.'6), the losses are nearly identical 
for both spectra in each case, and as we are con- 
cerned with flux ratio changes with wavelength, 
those losses do not affect our results. 



3 http: / /www. eso.org/gen-fac/pubs/astclim/lasilla/ 
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3. Data Analysis Methods and Uncertain- 
ties 

3.1. Continuum microlensing measure- 
ment 

The method we use to untangle microlens- 
ing and extinction is based on the measurement 
of the offsets between the continuum and the 
emiss i on line flux ratios (see, e.g 



19901 iMotta et all [200l IWucknitz et al 



Vandcrricst 



Wisotzki et al. 2003t iMediavilla et alJliooa 



2003; 



2011 



Sluse et al. I 120111 ). The multicomponent na- 
ture of quasar emission-lines imply that the 
emission-line spectra are produced over a wide 



(see e.g. Sulent 


ic et al. 


iMarziani et alJ ( 


2010h t 



2000h . According to 



(L1L) and the core of the high ionization lines 
(H1L) will be dominated by a component (FWHM 
~ 600 — 5000kms _1 ) that corresponds to the 
region of reverberation mapping typically large 
enough as to be insensitive to microlensing by so- 
lar mass objects. The broad wings of the emission 
lines (FWHM ~ 10000 km s _1 ), however, could 
arise from the inner parts of the BLR and may 
be microlensed. For this reason, we prefer to 
use exclusively the line cores (dominated by the 
NLR and the outer regions of the BLR) as ref- 
erence to set the baseline for no microlensing. 
To compute the core flux without attempting 
an an alytical decompositio n into several compo- 



nents ( Marziani et al. 2010h . we have used a nar- 



row band decomp osition similar to that used by 



Sluse et al. 



( 201 lh . Specifically, we define as core 
flux the continuum subtracted flux integrated in 
a relatively narrow velocity interval (from 25 to 
90 A depending on the line profile shape for the 
different sources) centered on the peak of the line. 
To accommodate the varying widths of the lines, 
the continuum estimate for each line requires win- 
dows with varying width as indicated below for 
each lens system. 

For each c omponent and each emission line we 
used DIPSO (|Howarthll2004h in STARLINK0 to fit 
a function y c = aX + b to the continuum on either 
side of the emission line, given a total wavelength 
range (Xa> X B ). The task also gives the error co- 



4 Support provided by the Starlink Project which is run by 
CCLRC on behalf of PPARC. 



efficients (Aa, Ab) in the continuum fitting. This 
error is largest in the bluest and reddest ends of 
the continuum, because it is affected by the re- 
sponse of the CCD. The flux under the continuum 
is then obtained as the integral below the fitted 
function y c , i.e. F c = (a/2)(X B ~X A ) 2 +b(X B -X A )- 
The error in the flux is estimated as AF C = 
(Aa/2){X B -X A f + Ab{X B - X A ). 

The emission line flux is obtained by integrat- 
ing the emission line profiles in each continuum- 
subtracted emission line using DIPSO. As com- 
mented above, we have separated the line core 
from the wings which could be affected by mi- 
crolensing. The error in the narrow emission line 
is estimated as the error in the continuum fit- 
ting. In those cases in which the emission line 
is affected by absorption lines, a narrower integra- 
tion window was chosen (10 to 15 A in the case of 
SDSS1029+2623). In most of the cases in which 
these absorptions are mild, they can be success- 
fully avoided. However, when the absorptions are 
broad and affect the central part of the emission 
lines (e.g. SDSS1029+2623) the measurements 
have correspondingly larger uncertainties. 

3.2. Impact of microlensing in the BEL 



.Nemirofl (1988) and ISchneider fc Wambsganss 
(|l990l ) suggested that, depending on the struc- 



ture of the BLR, mic rolensing coul d modi fy the 



broad line profiles. Abaias et al.1 ( 20021 ) have 



estimated which gravitational lens systems are 
more likely to show BLR changes due to mi- 
crolensing. Some examples of these variations 
i n the BLR have be e n presented bv |Filippenko 



1989); Chartas et al. (2002, 2004); Richards et al 



2004); iGomez- Alvarez et all (120061 ); ISluse et al 



20111 ) in the cases of MG0414+0534, H1413+117, 
SDSS1004+4112, and Q2237+0305. 

Superposition of the spectra for each image pair 
(see Figures 1, 6, 7, 11, 13, and 16) shows excel- 
lent matches between the emission lines profiles of 
HE1104-1805, SDSS1029+2623, Q0957+561 and 
HS0818+1227. Significant microlensing of BELs 
is detected only in the wings of one of the sys- 
tems, SDSS1004+4112 (see JO}. A slight en- 
hancement of the red wings that may be tenta- 
tively related to microlensing has been also de- 
tected in HE1104-1805 (see $L2]). Thus, the im- 
pact of microlensing on the BEL looks negligible 
except in SDSS1004+4112. The excellent matches 
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between the emission line profiles for each image 
pair imply that, except for fluctuations due to ab- 
sorptions or noise, the choice of the size of the 
line core has no impact on the results. To show 
this explicitly, we have compared the A-B mag- 
nitudes (averaged on all the lines for each sys- 
tem) computed using only the line core or the 
whole line (as estimated from the factor used to 
match the line profiles). We find differences of 
\{A — B) cores — (A — B) w hoie-ii7ies\—0-02, 0.03, 
0.01, and 0.06 for HE1104, SDSS1029, QSO0957, 
and HE0818, respectively. Thus, we find that for 
all the systems except SDSS1004-I-4112 the A-B 
emission line ratios do not significantly depend on 
the choice of the core width. In any case, we have 
used the core of the emission lines exclusively to 
estimate the emission line flux ratios (see section 

GEH). 

3.3. Estimate of accretion disk parameters 

For those cases in which chromatic microlens- 
ing is detected, we can study the structure of the 
accretion disk in the lensed quasar by estimating 
its size and temperature profile. We model the 
accretion disk as a Gaussian, / oc exp(— R 2 /2r 2 ), 
with radius variable with wavelength, r s cx X p . To 
estimate the probability of reproducing the mea- 
sured microlensing magnifications we have ran- 
domly placed a Gaussian source on microlens- 
ing magnification maps of 30 x 30 Einstein Radii 
squared (1000 x 1000 pixels) for SDSS1004+4112 
and 58.8x58.8 Einstein Radii squared (2000x2000 
pixels) for HE1104-1805 computed for each im- 
age using the Inverse Polygon Mapping method 
( Mediavilla et al. 2006). The convergence (k) and 
shear (7) for each image are select ed from available 



models in the literature (se e e.g. iKochanek et al 



2006t iMediavilla et al-lExih . We take a = 0.1 for 



the fraction of mass in compact objects, a reason- 
able value according to current e stimates (see e.g 




20021: IMediavilla et al 



Schechter & Wa mbs: 

20091 : IPoolev et all 1200911 We consider 1 M & mi 
crolenses. Following a Bayesian approach as in 
IMediavilla et~ai1 (|201lh . we estimate the proba- 
bility of r s and p conditioned on the measured 
microlensing magnifications for both uniform and 
logarithmic priors on r s . We have considered these 
two priors to analyze the sensitivity of our study to 



the treatment of the size pr ior (see iMorgan et al 



of 1 to 15 light-days (2.6 - 38.9 x 10 15 cm) for r s 
and a range of to 3 for p. In Sj4] we will ap- 
ply this method to SDSS1004+4112 and HE 11 04- 
1805. The results for r s and p are given with la 
errors. 

3.4. Dust extinction fitting 

Each lensed QSO image follows a different path 
through the lens galaxy, encountering different 
amounts of dus t and gas t h at pr oduce differen- 
tial extinction. Falco et all (|l999l ) measured the 
mean differential extinction, AE(B — V), in 23 
lens galaxies using HST broad-band filters. This 
extinction can affect not only the co ntinuum flux 



ratio but also the emission-line fluxes ( Motta et al 



2002t IMediavilla et~alll2005l 120091 |201l|). ThusT 
considering that the cores of the emission lines are 
affected neither by microlensing nor by intrinsic 
variability, measuring the emission lines flux ratio 
in several wavelengths provides us with a method 
to determine the existence of dust extinction in 
the system. We fitted the extinction curve to 
the magnitude difference in emi ssion lines for im- 
ages 1 and 2 using the equation ( Falco et al. 1999; 



Mufioz et al. 2004) 



2010l : Mediavilla et al-llioil ) . We consider a range 



mi(A)-m 2 (A) = -2.5 log +( J E 1 -.E 2 ) J Ry 

where M1/M2 is the constant magnification ra- 
tio, E\ — E% — AS is the extinction difference, 
and itV[A/(l + zl)] is the extinction curve in 
the lens rest frame. We minimized \ 2 P er num- 
ber of degree of freedom (xdof)- I* 1 the major- 
ity of the systems (except Q0957+561) we have 
only a few narrow emission lines in the optical 
par t of the spectra. We make our esti mates with 
the lCardelli. Clayton, fc Mathis (|l989l) extinction 
curve of the Milky Way (i.e. we fixed the param- 
eter Ry = 3.1) at the redshift of the lens galaxy. 
As is standard and to facilitate comparison of our 
results with those of other authors, magnitude dif- 
ferences are shown as a function of inverse wave- 
length in microns in the lens galaxy rest frame. 

3.5. Contamination from other sources of 
chromaticity 

There are two effects, intrinsic variability and 
contamination by the lens galaxy, that can pro- 
duce chromatic variations in the flux of lensed 
QSOs and, hence, mimic microlensing. 
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The continuum flux variation in QSOs is a well- 
know effect that does not significantly affect the 
NEL fluxes (|Petersonlll993l ). Intrinsic continuum 
variability combined with the time delay between 
images can produce a change in the flux ratios be- 
tween images that can be wavelength dependent, 
thus inducing changes in the chromaticity. These 
changes should be avoided if possible or at least 
estimated. In two of the objects, Q0957+561 and 
HE1104-1805 we can use data taken at two differ- 
ent epochs separated by the time delay to avoid 
the problem of intrinsic variability. 

In all the objects, we can e stimate the effects of 
intrin sic variability (following lYonehara. Hirashita fc 
20081 ) using the structure functio n inferred from 
the SPSS imaging data of quasars (jVanden Berk et al 
2004tllvezic et al.ll2004l ). We will consider the less 
favorable case; an intrinsic magnitude of Mj = 
—21 f or the quasar (jYonehara. Hirashita fc Richter 
2008), the bluest photometric band to measure 
variability and the two bands with the largest sep- 
aration in wavelength to estimate the chromaticity 
variation. In the case of SDSS1004+4112A,B with 
a measured time-delay of about 40 days, the ex- 
pected intrinsic variability is < 0.1 mag and the 
chromaticity change is < 0.03 mag. For HE1104- 
1805 with a measured time delay of about 150 days 
and for HS0818+1227 with a comparable theoret- 
ical delay, variability of ~ 0.1 mag and chromatic- 
ity change < 0.05 mag are predicted. Finally, for 
the largest separation systems, Q0957+561 and 
SDSS1029+2623, variability of < 0.2 mag and 
chromaticity change of < 0.08 mag are expected. 
Thus, changes in chromaticity, that are most sig- 
nificant to study the quasar structure, are rather 
small. 

The expected values of the intrinsic variabil- 
ity are in reasonable agreement with the analysis 



intrinsic variability in Q0957+561 and HE1104- 
1805. These effects are within the uncertainties for 
SDSS1004+4112 and, likely, for HS0818+1227 (al- 
though we lack on a measured time delay for this 
object). Finally, SDSS1029+2623 can potentially 
have relatively strong effects (~ 0.1 mag change 
in chromaticity) induced by intrinsic variability. 

On the other hand, in the cases of HS0818+1227 
and Q0957+561 the lens galaxy is bright and very 
close to one of the components on the sky (0'.'6 
and l'/O respectively) and some of the spectra 
may suffer contamination from the continuum of 

the lens galaxy. This continuum contamination is 

RichtetJ onger at longer wavelengths, but it does not af- 
fect the emission line fluxes. In these two cases, to 
avoid the continuum flux contamination we have 
considered the broad-band flux ratio obtained by 
CASTLE^I using HST imaging, in which the lens 
galaxy was modeled and subtracted. 



4. Results 

4.1. SDSS1004+4112 

SDSS1004+4112 is a five- image lens s ystem 
at z s — 1.734 discovered by llnada et"aL J2003I) 
with distances between components ranging from 
3'/7 to 14"6. The lens is a cluster a t zi = 0.68 

20081) . wh i ch ha s 



of lightcurves for Q0957+561 (IGoicoechea et al 
2008; IGoicoechea 1 120021 lOvaldsen et al.ll2003allbh 



SDSS1004+41 12 (IFohlmeister et al] 120081) and 
HE1104-1805 (IPoindexter et al.ll20071). In the case 
of HE1104-1805. IPoindexter et al l (|2007h specifi- 
cally studied the effect of intrinsic variability on 
flux ratios finding a global displacement of 0.1 
mag in magnitude differences due to the time de- 
lay without apparent changes in chromaticity in 
the optical (from the J to B photometric bands). 

In summary, we used photometry corrected for 
known time delays to avoid the effects induced by 



(|Oguri et al.1 12004 llnada et al 



Ota et all l2006h 



also been studied in X-rays 
This syste m has known CIV broad-line profile 
variations (jRichards et al.ll2004l) t hat are argued 



to ari s e either from microlensing (IRichards et al 



2004; iGomez- Alvarez et all 120061 : lAbaias et al 



20071 ) or due to small line-of-sight diff erences 
throu gh the quas ar absorbing outflows (iGreen 
20061 ). Recently, IFohlmeister et al.l (|2008l ) have 
measured a time delay of 40.6 ± 1.8 days for im- 
ages A and B, and 822 ± 2 days for C and D, 
detecting microlensing variability with an ampli- 
tude of the order of 0.1 5 mag between A and B 
(jFohlmeister et al.l 120081 ). 

Comparing the A and B spectra taken with 
the MMT we notice an enhancement in the blue 
wing and a decrement in the red wing of the CIV 
and SIV emission lines (Figure [I]). Lya and CIII] 
emission lines show smaller differences. Hence, our 



5 Cf A- Arizona Space Telescope LEns Survey, Kochanek, 
C.S., Falco , E.E., Impey, C, Lehar, J., McLeo d, B., Rix 
H.-W., http://www.cfa.harvard.edu/glensdata/ 
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results are consistent with Richards et al. I too4) 



Gomez- Alvarez et al 



(200G) 



and Lamer et al 



( 20061 ) although the amplitude of the enhancement 
of the blue wing is smaller than that observed pre- 
viously. This can be appreciated in Figure [5] where 
the A and B CIV emission line profiles taken in 
2004 with the Keck telescope are presented (data 
kindly provided by G.T. Richards). While the B 
component and the red part of the A component 
are basically the same in both epochs, the blue 
wing enhancement of component A is significantly 
smaller in 2008. This variability is the kind of 
gradual change in the line profile expected from 
microlensing. According to §3.1[ in what follows 
we will use the cores of the emission lines to com- 
pute flux ratios avoiding the effects of microlensing 
in the blue wings. 

A — B magnitude differences in the contin- 
uum and in the emission lines estimated from 
our spectra or obtained from the literature are 
shown in Figure [3] (see also Tables [2] and [3]). 
The A — B magnitude differences corresponding 
to the emission lines show no trend with wave- 
length (within uncertainties) and are distributed 
around < A — B >= —0.52 ± 0.07 mag sup- 
porting the absence of dust extinction and defin- 
ing the baseline for no microlensing magnifica- 
tion. In 2004 the continuum difference curve ob- 
tained from the spectra matched within errors 
the zero microlensing baseline defined from the 
low ionization emission lines. With small off- 
sets the broad-ban d ba sed continuum data from 
Oguri et all (|2004l ) and llnada et al.l (|2003l l2005h 
also match the baseline for no microlensing. This 
lack of microlensing evidence in the continuum in 
2004 (as the counterpart of the blue wing enhance- 
ments) was considered a serious drawback to in- 
terpret the enhancements in terms of microlensing 
(jGomez- Alvarez et al.l 12006). 

On the contrary, our A — B continuum differ- 
ence measurements (see Figure [3]) based on spec- 
tra taken in 2008, strongly depart from the zero 
microlensing baseline with an increasing trend to- 
wards the blue that wou ld include the X-r ay mea- 
surements obtained by lOta et all ([20061 ) . The 



magnitude difference in the continuum is consis- 
tent with CASTLES broad-band data. 

The A-B continuum differe nces corrected for 
the t ime delay measured by iFohlmeister et aL 
(|2008j ) change in the sequence: — 0.460±0.005 mag 



(2003-04), -0.283±0.007mag (2004-05), -0.339± 
0.005 mag (2005-06), and -0.381 ± 0.007 mag 
(2006-07). The lowest value, -0.46 ± 0.005 mag, 
is close to the mean magnitude difference in the 
emission lines, —0.52 ±0.07 mag, likely indicating 
that at this epoch (2003-04) the system showed 
little microlensing. 

Figure [4] shows a linear fit to the continuum 
data and the average of the emission line data. 
The magnitude difference variation in the contin- 
uum data (with a slope of 0.13 ± 0.04 mag //m -1 ) 
implies differences with respect to the emission 
lines of ~ 0.2 and ~ 0.5 mag at 7680 and 3320 
A respectively. Our results are consistent with 
the trend indicat ed by the X-ray continuum data 
( Ota et al. 20061 ). In summary, our data indi- 



cate negligible dust extinction and evidence of 
chromatic microlensing affecting the continuum. 
These results and the variability detected in the 
emission line profile give strong support to the hy- 
pothesis of microlensing to explain the enhance- 
ment in the blue wings. 

The structure of the accretion disk was studied 
using the procedure explained in section 13.31 We 
used the values (ka = 0.48, 7a = 0.59) and (kb = 
0.48, 7s = 0.48) taken from iMediavilla et al 



(2009) to obtain the magnification maps for the 
A and B images respectively. Applying this pro- 
cedure to microlensing measurements at three dif- 
ferent wavelengths corresponding to our MMT 
data (see Table 2]), we obtained the 2D prob- 
ability density functions (pdfs) shown in Fig- 
ure [5] for both linear and logarithmic grids in 
r s . From these distributions we obtain estimates 
r s = 7±31ight-days(18.1±7.8xl0 15 cm) andp = 
1.1 ± 0.4 for the linear prior and r s — 6^ 3 light — 



days(15.5 



+ 10.4 
-7.8 



x 10 i5 cm) and p = 1.0 ± 0.4 for 



the logarithmic prior. Although the value of p is 
consistent within uncertainties with the thin disk 
theory it is interesting to mention the trend in 
this and in other obje c ts to have p < 4 / 3 (se e 



Mediavilla et al.l l201ll: iBlackburne et al.l 1201 ll) 



The microlensing estimate for the size also exceeds 
substantially the estimate obtained from thin- 
disk theory (r,, ~ 0.31ightda ys = 0.78 x 10 15 cm, 



disk theory (r„ ~ U^gntda; 
Mosquera k Kochanekll201ll ) 



To study the impact of intrinsic variability in 
these results, we can compare the A-B difference 
we measured using the emission lines or the con- 
tinuum at 12500 A (see Table 4), where microlens- 
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ing and dust extinction should be less significant. 
We find a difference between both measurements 
(which is a conservative upper bound to contin- 
uum variability) of 0.08 mag, for an insignificant 
impact on the estimate of r s and p. 

4.2. HE1104-1805 

HE1104-1805 was discovered bv lWisotzki et al 



(|l993l ); it consists of two lensed images A and 
B separated by 3" 15 a t z R = 2.319. The lens 
galaxy was detected by ICourbin et al. ( 19981 ) at 
zi = 0.729. Image A is l'/l from the main lens 
galaxy. Variability in the continuum was de tected 



in spe ctra taken by IWisotzki et al.l (|1995) (op 



tical), Courbin et all ( 19981 ) (ne ar infrared), and 



Chartas et all (|2009f ) (X-ray). IPoindexter et al 



(|2007l ) monitored the system between 2003 and 
2006, concluding that the magnitude difference in 
the optical bands has changed from —1.7, when 
the lens was discovered, to —1.2 in their optical 
data (2006). These authors also provide a time 
delay estimation of 152. 2+^0 (1 o") days. 

The data obtained with the MMT and VLT 
show that the emission line profiles of both im- 
ages, A and B, are very similar, although some 
slight but interesting differences can be found in 
the broad components of CIV and SilV (Figures 
|6] and [7]) . Lya is only seen in our MMT spectra. 
The profile of the Mgll emission line is asymmetric 
both in A and B. CIII] presents heavy absorption 
lines both in the BEL and NEL. In the higher SNR 
data obtained from VLT it is clearly seen that the 
A spectrum shows several absorption lines, none 
of them present in the B spectrum. The profiles of 
CIV and SilV emission lines show a slight enhance- 
ment in the red wing of A compared to those of 
B both in MMT and VLT data. These wing en- 
hancements present only in high ionization lines 
might be evidence of microlensing. 

Figure [5] (see also Table [S]) presents the mag- 
nitude differences in the continuum and in the 
emission lines. We have also included data 
from the literature (Table [5]). The mean B-A 
magnitude difference (< B — A >= — 1.13 ± 
0.02 mag) corresponding to the emission lines 
obtained from MMT and VLT spectra is con- 



sisten t with the values deri ved by LW isotzki et al 



(jl995l) ( 1.14 mag) and ICourbin et al.l (l2000l) 

(— 1.16±0.04 mag). These values are also in agree- 
ment with the value estimated from infrared data 



-1.13 ± 0.03 mag (|Poindexter et alj|2007l) . These 
results confirm that the cores of the emission lines 
are not affected by microlensing and that little 
extinction is present. 

The B — A magnitude differences in the con- 
tinuum obtained from the MMT and VLT spec- 
tra show a slope that is in agreem ent with optical 



broad-band data obtained in 2006 ( Poindexter et al 



2007)F|. Broad-band data obtained several years 


before (Falco et al.ll999 


; Lehar et al.l2000l: Courbin et al. 


l2000t ISchechter et all 


2003^ are all consistent 



(slope — 0.16 ±0.03 mag /im : ) but are very differ- 
ent from our own recent data and that obtained 



bv IPoindexter et all (|2007l) . 

Linear fits to the magnitude differences of con- 
tinua are shown in Figure |H The slope for the 
magnitude differences in the emission lines and the 
IR data is 0.00 ± 0.06 mag ^m _1 which is consis- 
tent with no extinction and it is in good agreement 
wit h results obta i ned fr om near-infrared spectra 
by ICourbin et ail (|2000|) (A E < 0.01) and those 
found bv iFalco et al.l (Il999h (AE = 0.07 ± 0.1) 
using broadband data. The continuum data from 
the literature arc fitted in two separate sets: 1992- 
1994 data with a slope of —0.16 ± 0.03 mag /mi" 1 



and t he more recent data from IPoindexter et al 
(|2007h with a slope of 0.08 ± 0.06 mag /mi -1 . 
The slope of the linear fit to our 2008 contin- 
uum data (MMT+VLT), 0.12 ± 0.02 mag jum -1 , 
is thus in good agreement with the slope of 
0.08 ± 0.06 corr e spon ding to the 2006 data of 



Poindexter et al.l (|2007l) . but remarkably different 
from the value corresponding to 1992-1994 broad- 
band data —0.16 ± 0.03 mag /im -1 . Thus, mi- 
crolensing in HE1104-1805 has induced an extreme 
change in continuum slope that needs explanation. 
The chromaticity during the 1992-1994 epoch with 
an increasing amplitude towards the blue leads to 
an straightforward interpretation in terms of the 
magnification of the dominant component A. To 
explain the slope of the continuum corresponding 
to 2006-2008 epoch (under common assumptions 
about the unresolved source structure) we need 
to combine chromatic microlensing in both A and 
B. For instance, we can consider the combina- 



6 As we cannot correct our data for time delay, we have con- 
sidered both the time-d elay corrected and unco rrected op- 
tical data obtained by IPoindexter et at (2007) (magenta 
pentagons). Notice also that the lens galaxy continuum is 
very faint, so it cannot contaminate our spectra. 
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tion of two events of magnification in both A and 
B with a progressive increase in the strength of 
the B event from 1994 to 2008. This is only a 
qualitative example and simulations are needed to 
consistently reproduce microlensing chromaticity 
in each of the two epochs. 

Following the procedure described in section 
3.31 we have used the detected microlensing chro- 
maticity to study the structure of the accretion 
disk in HE1104-1805. We have done this for 
three sets of data: our VLT contin uum data 
from 2008, the IPoindexter et al.l (|2007t ) data cor- 



see iPoindexter et al.l 12008). 



rected for time delay and t he average of the dat a 
fro 



Falco et al. (1999), 



. - .Corn-bin et al.1 (119981), 
Lehar et al.l (|2000l) . and ISchechter et al.l (|2003l ) 
that consistently follow a common trend with 
wavelength. In Table [5] we present the microlens- 
ing measurements for each dataset. To compute 
the microlensing maps we have used the follow- 
ing projected densities and shears for each lens 
image (ha — 0.64, ^a = 0-52) and (kb = 0.3 3, 
7_b = 0.21) according to iMediavilla et al.l (|2009l ) 

The resulting pdfs are plotted in Figure IT0| 
and the expected values and uncertainties in Ta- 
ble [7] The pdfs corr esponding to the MMT /VLT 
(FigfTOji) and to the IPoindexter et al.l (|2008l ) data 



(FigflOb) are not as concentrated near the maxi- 
mum of the pdf as in the case of the broad-band 
data (Fig fTOb). These pdfs may present a sec- 
ondary maximum (perhaps due to the complex- 
ity of the microlensing phenomenon correspond- 
ing to this epoch) and, individually considered, 
are not very conclusive. However, the product pdf 
strongly increases the concentration of the prob- 
ability near the maximum and the significance 
of the estimates: r s — 6 ± 2 light — days(15.5 ± 
5.2 x 10 15 cm), p = 0.7 ± 0.1 for the linear prior 
and r s = 6±l light - days(15.5±^ x 10 15 cm), 
p = 0.7 ± 0.1 for the logarithmic prior. Our r s 
estimates correspond to one half light radius at 
the central wavelength of B filter, i? 1 / 2 (A4311) = 
8 ± 21ightdays(76 i 19 x 10 17 cm). This value 
is i n good agreement with the results obtained 
bv iMunoz et al.l ( 2011 ) with HST data and by 
Poindexter et al.l (|2008l ) from photometric moni- 
toring. The values of p are, however, considerably 
smaller. 

The microlensing-based size estimates are sig- 
nificantly larger than those inferred from the 
black-hole mass or from the observed I-band flux 



4.3. SDSS1029+2623 



SDSS1029+2623 was discovered bv llnada et al 



( 20061 ): it consists of two images A and B sep- 
arated by 22.5" at z s = 2.197 a nd a cluster 
lens g alaxies at Zl ~ 0.55. Recently lOguri et~al 



(2008) found a third image C 1.8" from B and 



several complex absorption systems in the emis- 
sion lines. 

Although the emission line profiles are similar 
for A and B (Figure fTTj) . there are several groups 
of absorption line systems affecting Lya and CIV 
that are associated with Mg l/Mgll/Fell absorp- 



tion systems, as found by lOguri et all (|2008l ) 



There are also self- absorption systems associated 
with Lya, SilV, and CIV lines that are present in 
both components but with significant differences. 
In spite of this we have attempted to determine 
flux ratios by defining suitable integration win- 
dows to avoid the absorptions. In the case of 
CIII], the emission line profiles are almost iden- 
tical in both components and show no absorption 
lines. Thus, the results derived from CIII] should 
be more reliable than the results inferred from the 
other lines. 

The B — A magnitude differences obtained from 
our data (continuum an d emission l i nes) com- 
pared to tho s e obta ined bv llnada et al. I (120061) and 
Oeuri et all |2008) are shown in Figure [TS] (Ta- 
ble |8|). O ur continuum flux ratio agrees well with 
the data cor responding t o the q and K broadband 
filters from lOguri et al.l (J2008) that were taken 
with the Keck at the same epoch. However, there 
is a difference of ~ 0.1 mag with the data taken 
at other epochs. This is explained by variability 
in the continuum between 2007 and 2008 (another 
peculi ar feature is that the measurement in the z 
band (Inada et al . 2006) is ~ 3 mag above all the 
other broad-band measurements) . In principle the 
variability could be attributed to microlensing or 
intrinsic variability of the quasar continuum com- 
bined with a time lag between both components. 
However, the strong chromaticity of the contin- 
uum flux ratio (of about 0.4 mag) that exceeds the 
~ 0.1 mag global offset between continuum flux 
ratios at different epochs excludes the explanation 
based on intrinsic continuum variability. Dust ex- 
tinction, on the other hand, cannot explain the 
chromaticity for the flux ratio inferred from radio 
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observations agrees with the flux ratio of the bluest 
continuum contrary to the expectations under this 
hypothesis. Thus, microlensing is the more likely 
explanation and is supported by the agreement of 
the flux ratios inferred from three of the lines, CIV, 
SilV, and Lya with the radio flux ratio. How- 
ever, the flux ratio inferred from the other emis- 
sion line, CIII] which presents the smoothest line 
profile, shows a large offset with respect to the 
baseline defined by the radio data that disagrees 
with the microlensing hypothesis. Under the hy- 
pothesis of chromatic microlensing we could follow 
the same steps as in the case of SDSS1004+4112 
to estimate the size and temperature profile of 
the quasar source in SDSS1029+2623. However, 
lens modeling in th is system is complex (see e.g. 
Kratzer et al.ll201ll ) and we defer this study to fu- 
ture work. 

4.4. Q0957+561 

The first known gravitational lens was dis cov- 
ered by IWalsh. Carswell. fc Wevmannl (|l979t k it 
has been studied in great detail. It consists of two 
images A and B with separation 6'.' 2. The source 
QSO is at z s = 1.41 and the main lens galaxy is at 
zi = 0.36 and is part of a poor cluster of galaxies. 
Comparison between A and B emission lines (Fig- 
ure [13]) do not show significant differences between 
the emission line profiles of CIV, CIII] and Mgll. 
This limits the possible impact of microlensing on 
the broad component of the emission lines. To 
quantify this impact we have compared the B/A 
flux ratios of the wings and the core of the CIV 
emission line (that has the highest S/N ratio) find- 
ing differences <10%. In any case we have com- 
puted flux ratios from the cores of the lines. 

Figure [H (see also Table EJ) shows the B - A 
magnitude differences in the continuum and in 
the emission lines for Q0957+561. This fig- 
ure also includes other data from the litera- 
ture and a re-a n alysis of HST/STIS data by 
Goicoechea et al.| ( 2005b ) . Ave raging the radio 

(|l992h at A = 6cm, 
at A = 13cm, and 
Haschick et al.l ()198lD at A = 6cm, the B - A 
magnitude difference uncontaminated by the lens 
galaxy continuum and free from dust extinction is 
obtained, < B - A > ra dio= 0.40 ± 0.03 mag. The 
B — A magnitude differences corresponding to the 
emission lines follow a decreasing trend towards 



data from 



_ .Conner et al 
Gorenstein et all (ll988F 



the blue compatible with extinction. A linear fit to 
the emission line magnitude differences (see Figure 
ITS"]) has a slope of — 0.04 mag //m -1 and a disper- 
sion of 0.09 mag. This dispersion is reasonable 
taking into account the intrinsic difficulty and the 
inhomogeneity of the data analysis procedures fol- 
lowed by the different authors especially regarding 
the criteria used to select the continuum. Towards 
the red this linear fit is fully consistent with the 
B — A radio measurements, confirming that the 
emission lines are not significantly affected by mi- 
crolensing. 

In Figure[l5]we also present an extinction curve 
fit to the B — A emission line magnitude dif- 
ferences (both from the literature and from our 
own measurements) . The best- fit parameters (ob- 
tained fixing the dust redshift to the lens red- 
shift) were: AE{B - V) = 0.02 ± 0.009 and 

R v = 2.0 ± 0.1 (xdof = L8 ) The data are 
also compatible with an extinction curve similar 
to the Milky Way (AE(B - V) = 0.02 ± 0.009, 
Xdof = 1-6)- It is remarkable that the fit- 
ting to the narrow emission lines is in agreement 
with other data obtain ed from the con tinuum 
AE(B-V) = 0.02±0.02 (|Falco et al.lll999l) . How- 
ever our results are not in agreement with the val- 
ues AE( B-V) = 0.068±0 005 and R v = 4.4±0.5 
found by iGoicoechea et ah ( 2005bl ). 

The continuum data also show a decreasing 
trend towards the blue with a slightly steeper 
slope (-0.05 ± 0.01 mag/mi" 1 for CASTLES and 
-0.05 ± 0.02 mag ^m" 1 for HST/STIS). In fact, 
the same extinction curve fitted to the emission 
line data with a shift of B — A ~ —0.4 mag fits 
well the continuum flux ratios from CASTLES and 
HST/STIS (notice that these continuum flux ra- 
tios obtained from HST data are not affected by 
lens galaxy contamination). This global shift be- 
tween the continuum and the baseline of no mi- 
crolensing magnification defined by the emission 
lines imply that the continuum is experimenting 
microlensing of ~ —0.4 mag and allows us to re- 
examine the microlensing history in Q0957+561 
based in this result. The differential B — A mag- 



nitude lightcurve of Q0957+561 (|Pelt et al.l 1 199a 



Oscoz et alJ 120 02) can be described as an event 



of 0.25 mag taking place from 1981 to 1986 and a 
quiet phase of mean value < B — A >^ —0.05 mag 
with fluctuations of less than 0.05 mag from 1987 
to 1999 (Oscoz et al. 2002). Previous attempts to 
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model the observed microlensing through simula- 
tions have not used the emission line flux ratios as 
microlensing zeropoint. These studies have either 
accepted any zeropoint for microlensing magnifica- 
tion, modeling microlensing variability of less than 
0.05 mag with respect to an un restricted zeropoint 
value (e.g. iRefsdal et al.ll2000f ) . or implicitly sup- 
posed that the zeropoint was placed at the mean 
value of the quiet phase (< B — A >~ —0.05 mag) 
modeling a microlensing amplitude o f less than 
0.05 mag (e.g. Wambsganss et all200Cll ). However, 
considering the zeropoint defined by the emission 
line flux ratios the correct procedure will be to 
model fluctuations of 0.05 mag with respect to 
a mean microlensing amplitude of < B — A >^ 
—0.4 mag. 

Although detailed microlensing simulations 
should be made to estimate physical parameters 
from the source and/or the microlenses, it seems 
that under this new perspective the likelihood of 
smaller sources (or bigger microdeflectors) will 
increase. 

On the other hand, there is an offset of ~ —0.2 
mag between the magnitude differences in the con- 
tinuum obtained with the MMT and those ob- 
tained with HST/STIS (data obtained at different 
epochs to correct for time delay). This difference 
between MMT and HST/STIS data might be ex- 
plained by (i) intrinsic variability (ii) a microlens- 
ing amplitude change, and (iii) continuum con- 
tamination by the lens galaxy. We plan to examine 
available photometric monitoring of Q0957+561 
covering 2008 (phased by the lag associated with 
component B) to ascertain the origin of this offset. 

If we compare the slopes of the linear fits cor- 
responding to emission lines and continua (CAS- 
TLES plus HST/STIS data) we found that the 
B — A offsets at Mgll and at OVI wavelengths 
have a difference of ~ 0.15 mag, these may be 
due to a wavelength dependence of microlensing 
(chromatic microlensing). This estimate is, how- 
ever, greatly affected by the uncertainties in the 
determination of the emission line flux ratios. 

In summary, for Q0957+561 our results indi- 
cate that: (i) there is no significant variation in 
the broad component of the emission line profiles, 
(ii) there is dust extinction affecting the emission 
lines and the continuum produced by dust likely at 
the same redshift as the lens galaxy and (iii) there 
is microlensing with amplitude B — A ~ —0.4 mag 



affecting the continuum. 

4.5. HS0818+1227 

HS0818+1227 was discovered byl agen fc Reimersl 
(2000); it consists of two lensed images A, B sep- 
arated by 2.6", with z s — 3.115 and z; = 0.39. 
Image B in this case is ~ 2 mag fainter than A. 
We scaled the continuum-subtracted B spectrum 
to match the emission line peaks in A (Figure 
[j"6]) . The A and B emission line profiles are very 
similar to each other and do not show significant 
differences in the BLR. 

Figure Q~7] shows the B — A magnitude differ- 
ences we calculated from the continuum (solid 
black squares) and from the cores of the emission- 
lines (solid black triangles) integrating our MMT 
spectra. The B — A magnitude differences cor- 
responding to the narrow emission lines (Table 
HOP will define a zero microlensing baseline of 
< tub — vriA >= 2.34 ± 0.03 mag. The average 
of CASTLES broadband data (< m B — rriA >= 
2.12 ± 0.03 mag) shows an offset of 0.22 mag 
with respect to this baseline. Note, however that 
the significance of this offset is dominated by the 
F555W data taken by CASTLES. According to 
section 3.5, part of this offset (0.1 mag) may arise 
from intrinsic variability. 

Ou r continuum data agree with Hagen fc Reimersl 
(|200dh . but do not match CASTLES (Tabled]), es- 
pecially in the reddest part. This discrepancy is 
due to the lens galaxy continuum. We estimate, 
using the integrated broad-band magnitudes ob- 
tained by CASTLES for the lens galaxy, that the 
contamination is ~ 40 % of the lens galaxy flux. 

Considering the emission lines, our results in- 
dicate negligible dust extinction and posible evi- 
dence of microlensing. 

5. Conclusions 

The method we use in this paper allows us to 
separate microlensing from dust extinction with- 
out a theoretical model for the lens system. We 
have demonstrated the method for the most com- 
plicated cases: doubly- imaged quasars. 

We tested the hypothesis that the cores of the 
emission lines do not vary with time by compar- 
ing our own magnitude differences in the emission 
lines with values from the literature that were ob- 
tained at different epochs, including values cor- 
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rected for measured time-delays, and we conclude 
that they are nearly constant in time. Thus, ex- 
cept in cases where extinction is significant, the 
magnitude differences in the emission line cores 
are reliable estimators of the intrinsic magnitude 
differences unaffected by microlensing. 

Following Yonehara. Hirashita fc Richterl ( 2008f) 
we have estimated the impact of time delays 
in our microlensing measurements for our ob- 
jects. In the worst case scenario (Q0957+561, 
SDSS1004+4142, and SDSS1029+2623) a time 
delay can introduce variabilities < 0.2 mag and 
chromaticities < 0.1 mag. The measurements we 
obtain for those objects are at least twice the esti- 
mated values. Although more data are needed for 
confirmation, it appears that time-delay induced 
variability has a modest impact. 

Differences in the wings of the CIV and SilV 
broad emission line profiles are found in A and B 
im ages of SDSS10 4+41 12, as detected previously 
bv lRichards et al~ ( 2004 ). but the enhancement in 
the blue wing is smaller than observed in 2004. 
In HE1104-1805 we also have detected a slight en- 
hancement in the red wings of CIV and SilV in 
image A with respect to image B. 

The average microlensing magnification free 
from extinction was obtained as the difference in 
magnitudes between the emission lines and the 
continuum. The latter was obtained directly from 
the spectra in those cases where there is no con- 
tamination by the lens galaxy (all systems ex- 
cept HS0818+1227 and Q0957+561), otherwise we 
used HST continuum data free from lens galaxy 
contamination available in the literature. Sig- 
nificant chromatic microlensing was detected in 
SDSS1004+4112, SDSS1029+2623, and HE1104- 
1805. 

Below is a summary of the results for each sys- 
tem: 

1. We detected a blue wing enhancement in 
the high ionization lines of SDSS1004+4112 
that are qualitatively similar to the ef- 
fect described by previous authors but of 
smaller amplitude. We have also detected 
strong chromatic variability in the contin- 
uum. The presence of variability in both 
lines and continuum supports the hypothe- 
sis of microlensing to explain the blue wing 
enhancements in the lines. Our data indi- 



cate negligible dust extinction. We infer an 
accretion disk size of r s = 7 ± 3 lightdays = 
18.1 ± 7.8 x 10 15 cm at A resi = 3363A and 
a wavelength dependence of the size with 
exponent p = 1.1 ± 0.4. 

In HE1104-1805 we find no extinction but we 
detected chromatic microlensing with large 
variations between two epochs that change 
the sign of the continuum slope. We esti- 
mate r s = 6 ± 2 lightdays = 15.5 ± 5.2 x 
10 15 cm at X rest = 3363A and p = 0.7 ± 0.1. 
This size is greater than those inferred from 
the thin disk theory and either the intrin- 
sic flux or the c e ntral black hole mass by 
Poindexter et al.l (2008). However, our re- 



sults agree with re cent values presented in 



Munoz et al.l(|2011l) rJA3363) = 7±4 light- 
days and p = 1.1 ± 0.6 The value of p we 
obtain is significantly smaller than the 4/3 
value predicted by the standard model. Sim- 
ilar discrepancies have been found in other 



objects (Poindexter et al. 2008; Flovd et al 



2009 ; Morgan et al.l 1201 0; Blackbur ne et al 



2011 



20111 ). 



Mediavilla et al.l |201 It I Munoz et al 



3. SDSS1029+2623 is affected by strong chro- 
maticity of about 0.4 mag that can be ex- 
plained neither by dust extinction nor by in- 
trinsic variability. Chromatic microlensing 
is the most probable explanation although 
not all the available data fit well within this 
hypothesis. 

4. A ICardelli. Clavton. fc Matins! (Il989h ex- 
tinction law with AE(B - V) = 0.02 ± 0.09 
and R v = 2.0 ± 0.1 (xdof = L8 ) can be 
used to fit, within uncertainties, both the 
continuum and the emission line flux ra- 
tios in Q0957+561. There is a global offset 
between the continuum and emission line 
B — A magnitude differences that implies 
a microlensing amplitude of < B — A >^ 
—0.4 mag. There are marginal indications 
of chromatic microlensing. 

5. We detect evidence of microlensing but 
not extinction (within uncertainties) in 
HS0818+1227. 
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Fig. 1.— Lya, SilV, CIV, CIII] emission line profiles for SDSS1004+4112 vs. observed A. The red line 
represents the continuum subtracted emission lines for A. The black line represents the continuum subtracted 
emission line for B multiplied by a factor to match the peak of A. The factors are shown in each panel. 
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Fig. 2. — CIV emission line profile comparison for SDSS1004+4112. Red and black lines rep resent A and B 
MMT spectra respectively, magenta and blue represent A and B Keck spectra obtained by iRichards et al 
(|2004l ) respectively. 
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SDSS1004+4112 




_Q Q I I I I I I I I L_U I I ILJ I I I ll I I_J I LU I I L_ 

1 2 3 4 5 6 

Fig. 3. — Magnitude differences niA —m-B vs Aq 1 (A in the lens galaxy restframe) for SDSS1004+4112. We use the standard 
units of /an -1 for extinction studies, which are convenient to cover the range of observ ed A. Sol id pentagons repre sent the 
integrated continu um obtained fro m (broad-band) CASTLE S {red), llnada et al,l d2003f) {green), llnada et al.l (I2005T) {blue), 
lOguri et alj d2004T) (magent a), and iFohlmcist er et al. | ^2008l> {cyan). The green open pentagon represents the X-ray data 
obtained bv lOta et al . (2006) (for display convenience, we shifted it in wavelength from 60 to 290 fan; i.e from 28 to 5.8 fim 
in the rest frame). The black and blue sguares represent the magnitude differences from the integrated continuum in our spectra 
{solid) and from the integrated fitted continuum under the emission lines {open) for two different exposures. Black and blue 
triangles are the magnitude difference in emission line core. 
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Fig. 4. — Model fitted to the data shown in Figure [3] Squares and triangles represent continuum and NEL 
data respectively. Black lines represent the function fitted to the continua and the average of the emission 
line cores. Dashed lines are the standard deviations for the continuum fits and the standard error of the 
mean for the emission line cores. 
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P P 

Fig. 5. — Two-dimensional pdfs obtained using the measured chromatic microlensing for SDSS1004+4112 
(Table SJ for both linear (left) and logarithmic (right) grids in r s . Contours correspond to 0.5er, la, 1.5a, 
and 2<T respectively. We estimate r s = 7 ± 3 light — days(18.1 ± 7.8 x 10 15 cm) and p — 1.1 ± 0.4 for the linear 
prior and r s = 6^3 light — days(15.5t7° 8 4 x 10 15 cm) and p = 1.0 ± 0.4 for the logarithmic prior. The dashed 
line corresponds to the value predicted by the thin disk model (p = 4/3) 
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Fig. 6.— CIV, SilV, CIII] emission line profiles for HE1104-1805 vs. observed A. Upper panel MMT 
spectra. The red line represents the continuum-subtracted emission lines for A. The black line represents 
the continuum subtracted emission line for B multiplied by a factor to match the peak of A. The factors 
are shown in each panel. Bottom panel same as upper panel but for VLT spectra. 
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HE1104-1805 



HE1104-1B05 




Fig. 7. — Lya, Mgll emission line profiles for HE1104-1805 vs. observed A. The red line represents the 
continuum subtracted emission lines for A. Black line represents the continuum subtracted Lyct emission 
line for B multiplied by 2.7 to match the peak of A (red line). 
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Fig. 8. — Magnitude differences rriA — ttib vs A^ 1 (A in the lens galaxy restframe) for HE1104-1805. Black and blue represent 
the magnitude differences obtained from MMT and VLT spectra respectively. Solid squares are the magnitude differences 
in the continuum, open squares in the integrated continuum under the emission line, and solid trianqles in the emission line 
core. The broadband da ta obtai ned from other autho rs are pl otted as pentagons in different colors repr e sentin g: CASTLES 
(rafl.lLehar et aI.H2000l) ( green), ICourbin et all J199ST) ( cyan) . ISchechter et al] d2003fl (fefae). lFalco et al.1 il999T) (open black), 
Poindexter_etjil. J2003) Spitzer IRAC (solid black). Magenta pentagons represent the optical broadband data obtained by 
Poindexter et aU (2007) with (solid) and without (open) time-delay correction. 
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Fig. 9. — Model fitted to the data shown in Figure [8] Black lines represent the fitted function to the continua 
(squares) for MMT and VLT data and the average of the emission line cores (triangles) respectively. The 
red line represents the fitted function to the broadband dat a in the literature (CASTLES; iLehar et al. 2000; 
ICourbin et al. 19981 Falco et al. 19991 : Schechter et al.112003 ) at the same epoch . The magenta line represents 
the fitted function to the broadband data obtained bvlPoindexter et al.l d2007l) with and without time-delay 
correction. The data obtained in the infrared by IPoindexter et al.l ( 2007 ) are plotted as black triangles. 
Dashed lines are the standard deviation for each fit (continua and broadband data) and the error of the 
mean for the emission line cores. 
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Fig. 10. — Two-dimensional pdfs obtained using the measured chromatic microlensing for HE1104-1805 (Table [6} for both 
linear (top) and logarithmic (bottom) grids in r s . Contours are 0.5<r, ltr, and l.5o~ con fidence levels respective ly. From left to 
right pdfs for: our MMT/VLT data (a), average of broadband data previous to 2003 (bl. lPoindexter et all {20071 data corrected 
by time delay (c), and the intersection among the three previous maps (d). In the intersection maps we also show the contour 
corresponding to 2cr confidence level. 
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Fig. 11— Lya, SilV, CIV, CIII emission line profiles for SDSS1029+2623 vs. observer A. The red line 
represents the continuum subtracted emission lines for A. The black line represents the continuum subtracted 
emission line for B multiplied by a factor to match the peak of A. The factors are shown in each panel. Sky 
lines are seen on both sides of CIII] . 
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SDSS1029+2623 




Fig. 12. — Magnitude differences mj — tua vs A 1 (A in the lens galaxy re stframe) for SDSS1029+2623. 
Solid squares repr esent the integrated continua, in color those obtained by llnada et all ((2006) (red) and 
Oguri et al. |2008) (green and magenta represent data obtained in 2007 and in 2008 respectively), and in 
black those obtained from our spectra. Open black squares represent the difference in the integrated fitted 
continua under the emission lines. Black triangles are the magnitude differences in the emission line cores. 
The blue line repres ents the magnitude difference and its error (blue dashed lines) at radio wavelengths 
(jKratzer et al.ll201ll ). 
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Q0957+561 Q0957+561 Q0957+561 




3500 3600 3700 3800 390 4400 4600 4800 6200 6400 6600 6800 700 

A (A) A (A) A (A) 

Fig. 13. — CIV, CIII], Mgll emission line profiles for Q0957+561 vs observed A. The red line represents 
the continuum subtracted emission lines for A. The black line represents the continuum-subtracted emission 
lines for B multiplied by a factor to match the peak of A. The factors are shown in each panel. 
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Fig. 14. — Magnitude differences ms — rriA vs A^ 1 (A in the lens galaxy restframe) for Q0957+561. Black and blue represent 
the magnitude difference obtained from MMT and HST spectra respectively. Solid squares are the magnitude difference in 
the integrated continuum, open squares are the integrated continuum under the emission lines, and solid triangles are the 
integrated emission line cores. The data obtained from other authors are also plotted following the previous code: solid squares 
are broadband data or integrated continuum, and solid tria ngles are em i ssion line cores. The code for the colors is: red 
CASTLES data for the continuum llBernstein et al1ll997t) and Vanderrics t| ljl993T l data f or Mgll emission line, cyan estimated 
from the spectra of Mediavilla et al. (2000), maoerato lGoicoechea et a.1.1 (l2005bTi . and greera lSchild fc Smithl l|l99ll) . The Red line 
is the mean magnitude difference and standard d eviation of the mean (red dashed lines) at radio wavelengths l|Conner et al] 
l!992l : lGorenstein et al. 1 119881 : lHaschick et~aLlll98lH . 
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Fig. 15. — Model fitted to the data shown in Figure 1141 The black lines represent the function fitted to the emission 
line cores (solid triangles) and the continua (solid squares) obtained with HST respectively. Dashed lines are the standard 
deviation for each fit. The blue curve represent the dust extinction function fitted to the emission line cores using Ry =3.1 
(E(B — V) = 0.02 ± 0.09 with Xdof = 1-6)- The green line is the dust extinction fitted using variable Ry (Ry = 2.0 ± 0.1, 
E(B — V) = 0.02 ± 0.09 with X 2 dof = 1-8)- The curves shifted —0.4 mag (dashed blue and dashed green) fit the HST continua. 
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HS0818+1227 HS0818+1227 HS0818+1227 
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A(A) A(A) X(k) 

Fig. 16. — OVI, Lya, and CIV emission line profiles for HS0818+1227 vs. observed A. The red line represents 
the continuum-subtracted emission lines for image A. The black lines represent the continuum subtracted 
emission lines for B multiplied by factors to match the peak of A. The factors are shown in each panel. 
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HS0818+1227 




Fig. 17. — Magnitude differences mj- w,a vs A 1 (A in the lens galaxy rest frame) for HS0818+1227 . 
Solid squares represent the integrated continuum obtained by CASTLES {red), lHagen fc Reimersl (|200nh 
( magenta), and from our sp ectra (black). The magenta triangle is the magnitude difference in CIV obtained 
bv lHagen fc Reimersl (2000). Open black squares represent the integrated fitted continua under the emission 
lines. Black solid triangles are the magnitude difference in the emission line cores. 
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Table 1 
Log of observations 



Objects 


Pair" 


A b (") 


Instrument 


Grating 


Date 


Airmass 


P.A. C 


Sccing d 


Exposure 6 


HS0818+1227 


AB 


2.6 


MMT/Bluc-Channcl 


300 


2008/01/12 


1.168 


-36.10 


0.59 


1800 


Q0957+561 


AB 


6.2 


MMT/Blue-Channcl 


300 


2008/01/12 


1.096 


168.51 


0.61 


900 




A 




HST/STIS 


G230L 


1999/04/15 








1900 




B 




HST/STIS 


G230L 


2000/06/02 








1900 




A 




HST/STIS 


G430L 


1999/04/15 








900 




B 




HST/STIS 


G430L 


2000/06/03 








900 




A 




HST/STIS 


G750L 


1999/04/15 








660 




B 




HST/STIS 


G750L 


2000/06/03 








660 


SBSS1004+4112 


AB 


3.8 


MMT/Blue-Channel 


300 


2008/01/12 


1.028 


200.40 


0.61 


2 x 900 


SBSS1029+2623 


AB 


22.6 


MMT/Blue-Channel 


300 


2008/01/11 


1.072 


11.12 


0.67 


1800 


HE1104-1805 


AB 


3.2 


MMT/Blue-Channel 


300 


2008/01/11 


1.766 


114.66 


0.67 


1000 




AB 


3.2 


VLT/FORS2 


300V 


2008/04/07 


1.315 


64.18 


0.60 


3 x 250 



a Pair or image observed 
b Scparation between images in arcsec 
c Position angle in degrees E of N 
d Sccing in arcsec 
°Scconds of time 



34 



Table 2 

Summary of Known Quasar Image Properties 



Lens Name z L a z s b Filter 1/A d (pm 1 ) Am (mag) e Source f 



HS0818+1227 



0.39 1.3115 



Q0957+561 



0.36 



1.41 



SBSS1004+4112 0.68 1.734 



F160W 


0.65 


2.17 ±0.03 


1 


F814W 


1.23 


2.11 ±0.00 


1 


R 


1.55 


1.2 


2 


(CIV) 


1.57 


2.26 


2 


F555W 


1.80 


2.07 ±0.22 


1 


radio 


0.00 


0.32 ±0.02 


3 


radio 


0.00 


0.4 ±0.2 


22 


radio 


0.00 


0.48 ±0.03 


23 


F160W 


0.65 


0.08 ±0.01 


1 


F814W 


1.23 


0.02 ±0.02 


1,4 


(Mgll) 


1.49 


0.48 ±0.04 


5 


(Mgll) 


1.49 


0.31 ±0.02 


6,7 


(Mgll) 


1.49 


0.22 


8 


R 


1.55 


-0.022 ±0.013 


9 


F555W 


1.80 


0.01 ±0.02 


1,4 


V 


1.83 


-0.077 ±0.023 


9 


(CIII]) 


2.17 


0.27 ±0.03 


5 


(CIV) 


2.68 


0.28 ±0.02 


5 


(NV) 


3.34 


0.43 ±0.07 


5 


(Lya) 


3.40 


0.15 ±0.03 


5 


F284M 


3.52 


0.06 ±0.04 


10 


F277M 


3.61 


0.02 ±0.03 


10 


F248M 


4.03 


0.06 ±0.04 


10 


F140LP 


4.55 


0.00 ±0.06 


10 


F160W 


0.65 


-0.47 ±0.04 


1 


z 


1.10 


-0.45 ±0.09 


11 


z 


1.10 


-0.45 ±0.06 


12 


F814W 


1.23 


-0.31 ±0.15 


1 


F814W 


1.23 


-0.34 ±0.12 


13 


i 


1.30 


-0.40 ±0.06 


11 


i 


1.30 


-0.40 ±0.06 


12 


r 


1.60 


-0.460 ±0.005 


14 


r 


1.60 


-0.283 ±0.007 


14 


r 


1.60 


-0.339 ±0.005 


14 


r 


1.60 


-0.381 ±0.007 


14 


r 


1.60 


-0.39 ±0.08 


11 


r 


1.60 


-0.39 ±0.06 


12 


F555W 


1.80 


0.09 ±0.17 


1 
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Table 2 — Continued 



Lens Name z^ a 


Z5 b 


Filter c 


1/A d (pm- 1 ) 


Am (mag) e 


Source f 






g 


2.08 


-0.37 ±0.08 


11 






g 


2.08 


-0.37 ±0.06 


12 






u 


2.84 


-0.40 ±0.08 


12 


SBSS1029+2623 0.55 


2.197 


radio 


0.0 


-0.272 ±0.013 


24 






z 


1.10 


0.44 ± 0.09 


15 






I 


1.24 


-0.06 ±0.02 


16 






i 


1.30 


0.02 ±0.04 


15 






R 


1.52 


-0.04 ±0.02 


16 






R 


1.52 


0.12 ±0.02 


16 






r 


1.60 


0.03 ±0.01 


15 






V 


1.82 


-0.05 ±0.02 


16 






g 


2.08 


-0.01 ±0.01 


15 






g 


2.08 


0.09 ±0.02 


16 






B 


2.25 


-0.17±0.02 


16 






u 


2.84 


-0.06 ±0.04 


15 


HE1104-1805 0.73 


2.32 


IRAC 8.0//m 


0.13 


-1.15 ± 0.01 


17 






IRAC 8.0fim 


0.13 


-1.11 ±0.01 


17 






IRAC 5.8fiia 


0.17 


-1.13 ±0.02 


17 






IRAC 5.8Atm 


0.17 


-1.09 ±0.02 


17 






IRAC 4.5/zm 


0.22 


-1.37 ±0.04 


17 






IRAC A.bfiia 


0.22 


-1.13 ±0.03 


17 






IRAC 3.6,um 


0.28 


-1.44 ±0.02 


17 






IRAC 3.6/im 


0.28 


-1.15 ±0.02 


17 






K 


0.45 


-1.35 ±0.11 


18 






"PI -1 r* ATI T 

F160W 


0.65 


1 A A I A AO 

— 1.44 ± 0.03 


1 






F160W 


0.65 


-1.47 ±0.03 


19,20 






J 


0.77 


-1.53 ±0.08 


18 






J 


0.80 


-1.45 


17 h 






J 


0.80 


-1.38 


17 1 






F814W 


1.23 


-1.55 ±0.03 


1 






F814W 


1.23 


-1.61 ±0.02 


20 






F814W 


1.23 


-1.63 ±0.06 


19 






I 


1.27 


-1.25 


17 1 






R 


1.55 


-1.35 


17 h 






R 


1.55 


-1.23 


17 1 






F555W 


1.80 


-1.78 ±0.03 


1 






F555W 


1.80 


-1.76 ±0.03 


20 






F555W 


1.80 


-1.82 ±0.05 


19 



36 



Table 2 — Continued 



Lens Name z^ a z$ h Filter c 


1/A d (^m- 1 ) 


Am (mag) c 


Source f 


V 


1.83 


-1.748 ±0.03 


21 


B 


2.28 


-1.23 


17 h 


B 


2.28 


-1.13 


17 1 



a Lens galaxy redshift 
b Lensed quasar redshift 

c Filter or, when available, the line emission flux between parentheses 

d Inverse of the central wavelength (rest frame) Radio wavelengths are approx- 
imated as fi m -1 in our plots 

e Magnitude difference of B-A. Except for SDSS1004+4112 and HE1104-1805 
where we show A-B 

h Magnitude difference with time delay correction 

'Magnitude difference without time delay correction 

" REFERENCES: (1) CASTLES; ( 2) iHaeen fc Reimersl (l2000f): (3 ) 



Conner et all (Il992h : (4) iBernstein et all (fl997l ); (5) iGoicoechea et all (|2005bl ); 
(6) ISchild fc Smith! (Il99ll ): (7) IVanderriestl (Il993l): (8) Mediavilla (privat e 
comunication based on data taken in 1 997); (9) iGoicoechea et al.l (2005al ): 



(10) iDolan et al.l dl995h : (11) llnada et all (120031) : (12) lOeuri et all (l2004h . (13) 
Inada et all (l2005h : (14) iFohlmeister et al.l d 2008h: CIS) llnada et all (12009): (16 ) 
Qguri et all (|2008|) ; (17) IPoindexter et al.l 120071) ; (18 ) ICourbin et all (1998 1) : 



(19) iLehar et alJ (l2000h: ( 201 iFalco et al 
(22) [Haschick et al.l 
( 201lh . 



( 211 ISchechter et alJ (12 003) ; 



(jl98lh : (23) iGorenstein et all (|l988l ): (24) iKratzer et all 
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Table 3 

sdss1004+4112 magnitude differences 



Region 


A c (A) 


Window (A) 


rriA — rtiB a (mag) 


rriA — ti%b (mag) 


Continuum 


3320 


3100-3700 


0.07 ±0.05 


0.03 ± 0.08 




3820 


3600-4040 


0.05 ±0.04 


0.01 ±0.05 




4230 


3970-4450 


-0.05 ±0.03 


-0.12 ±0.03 




4500 


4350-4750 


-0.04 ±0.04 


-0.08 ±0.04 




5215 


4600-5550 


-0.01 ±0.03 


-0.09 ±0.04 




6650 


6400-6870 


-0.27 ±0.03 


-0.28 ±0.04 




7680 


7150-8100 


-0.25 ±0.05 


-0.25 ±0.05 


Line 


Lyal216 


3310-3345 


-0.45 ±0.05 


-0.54 ±0.08 




S1IVA1400 


3790-3870 


-0.56 ±0.04 


-0.57 ±0.05 




CIVA1549 


4220-4260 


-0.57 ±0.03 


-0.61 ±0.03 




HeIIA1640 


4470-4510 


-0.66 ±0.04 


-0.53 ±0.04 




CIII]A1909 


5210-5250 


-0.38 ±0.03 


-0.41 ±0.04 




CIIA2326 


6630-6680 


-0.52 ±0.03 


-0.48 ±0.04 




MgIIA2800 


7640-7760 


-0.30 ±0.05 


-0.33 ±0.05 



a Exposure 1 
b Exposure 2 
integration window. 
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Table 4 

sdss1004+4112 chromatic microlensing 



A c (A) 


Amc — Am_L a (mag) 


3700 


0.60 ±0.02 


6338 


0.40 ±0.02 


12500 


0.08 ±0.04 



a Diffcrcncc between the mag- 
nitude difference in the contin- 
uum and in the emission fines 
(m B - m A ) c - (m B - m A ) L 
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Table 5 

he1104-1805 magnitude differences 



Region 


A c (A) 


Window (A) 


m A -m B a (mag) 


m A - m B h (mag) 


Continuum 


4037 


3800-4350 


-1.10 ±0.04 






4638 


4400-4900 


-1.13 ±0.01 


-1.11 ±0.02 




5143 


4850-5400 


-1.22 ±0.01 


-1.16 ±0.02 




6338 


5600-6800 


-1.15 ±0.02 


-1.24 ±0.01 




9293 


8700-9600 




-1.32 ±0.03 


Line 


Lycd216 


4013-4050 


-1.02 ±0.04 






S1IVA1400 


4600-4670 


-1.12 ±0.01 


-1.12 ±0.02 




CIVA1549 


5080-5160 


-1.22 ±0.01 


-1.18 ±0.02 




CIII]A1909 


6250-6360 


-1.20 ±0.02 


-1.04 ±0.01 




MgIIA2800 


9240-9290 




-1.21 ±0.03 



a MMT data 

b averaged VLT data 

integration window. 
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Table 6 

he1 104-1805 chromatic microlensing 



A c (A) 


Amc — AniL a (mag) 


4380 


-0.09 ±0.02 


6470 


-0.21 ±0.02 


12500 


-0.34 ±0.04 


5550 


0.65 ±0.03 


8140 


0.46 ±0.03 


15500 


0.32 ±0.03 


3700 


-0.15 ±0.02 


6338 


0.16 ±0.02 


12500 


0.31 ±0.02 



a Difference between the mag- 
nitude difference in the contin- 
uum and in the emission line core 
(niB~mA)c-{mB-mA)L- Mea- 
surements corresponding to our 
MMT/VLT data, average data 
from before 2003 (see text), and 
IPoindexter et al. (2007) data cor- 
rected by time delay respectively. 
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Table 7 

he1 104-1805 accretion disk parameters (1(7 error) 







Linear 


Logarithmic 


Data 


P 


r s x 10 15 (cm) 


P 


r s x 10 15 (cm) 


This work 


1.7 ±0.8 


28.5 ± 10.4 


1.7 ±0.8 


23.3 ±5.2 


Average Lit. 


1.1 ±0.7 


23.3 ± 10.4 


0.9 ±0.2 


10.4 ±5.2 


Poindexter et al. (2007) 


0.8 ±0.3 


23.3 ± 7.8 


1.1 ±0.6 


18.1 ±5.2 


Intersection 


0.6 ±0.1 


15.5 ±5.2 


0.7 ±0.1 


15.5 ±2.6 
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Table 8 

SDSS1029+2623 magnitude differences 



Region 


A c (A) 


Window a (A) 


m B - m A (mag) 


Continuum 


3888 


3650-4450 


-0.07 ±0.04 




4466 


4250-4670 


-0.01 ±0.02 




4952 


4700-5150 


0.05 ±0.02 




6103 


5750-6280 


0.10 ±0.01 


Line 


Lyal216 
SHVA1400 
CIVA1549 
CIII]A1909 


3890-3905 
4450-4470 
4965-4975 
6090-6115 


-0.39 ±0.04 
-0.40 ±0.02 
-0.33 ±0.02 
-0.03 ±0.02 



integration window. 
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Table 9 

q0957+561 magnitude differences 



Region 


A c (A) 


Window c (A) 


m B - m A a (mag) m B - m A h (mag) 


Continuum 


2487 


2400-2580 




-0.15 ±0.05 




2930 


2750-3100 




-0.20 ±0.06 




2990 


2750-3100 




-0.20 ±0.06 




3370 


3210-3600 


-0.39 ±0.02 


-0.12 ±0.04 




3730 


3500-3890 


-0.35 ±0.01 


-0.14 ±0.03 




3950 


3850-4130 


-0.33 ±0.01 


-0.13 ±0.03 




4600 


4270-4850 


-0.30 ±0.006 


-0.10 ±0.01 




6760 


6230-7000 


-0.23 ±0.01 


-0.01 ±0.03 


Line 


OVIA1032 


2490-2515 




0.32 ±0.05 




Lyal216 


2920-2945 




0.13 ±0.06 




NVal240 


2980-3000 




0.30 ±0.06 




SHVA1400 


3355-3400 


0.28 ±0.02 


0.17 ±0.04 




CIVA1549 


3710-3755 


0.38 ±0.01 


0.40 ± 0.03 




HeA1640 


3935-3965 


0.18 ±0.01 


0.20 ±0.03 




CIII]A1909 


4560-4630 


0.42 ± 0.006 


0.19 ±0.01 




MgIIA2800 


6730-6785 


0.30 ±0.01 


0.26 ±0.03 



a MMT data 
b HST data 
integration window. 
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Table 10 

HS0818+1227 magnitude differences 



Region 


A c (A) 


Window a (A) 


m B - m A (mag) 


Continuum 


4245 


3800-4800 


1.96 ±0.07 




5004 


4700-5070 


2.0 ±0.1 




6374 


6000-6650 


1.36 ±0.02 


Line 


OVIA1037 


4220-4270 


2.35 ±0.07 




Lyal216 


4985-5015 


2.4 ±0.1 




CIVA1549 


6340-6390 


2.40 ±0.02 



a Integration window. 
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